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We calculate analytically the low temperature quasi-particle scattering rate, the conductivity, 
and the specific heat in weakly disordered metals close to a quantum critical point, via the use of 
a proper fluctuation potential V{q, oj) between the quasi-particles. We obtain typical Fermi liquid 
results proportional to and T respectively, with prefactors which diverge as power laws of the 
control parameter a upon approaching the critical point. The Kadowaki-Woods ratio is shown to 
be independent of a (possibly times a logarithmic dependence on a) only for the case of three- 
dimensional ferromagnetic fluctuations. Our results are consistent with experiments on the eight 
materials CeCoIns, Sr 3 Ru 207 , YbRh 2 Si 2 , La 2 -a:Cea:Cu 04 , Tl2Ba2Cu06+2;, CeAuSb 2 , YbAlB4, and 
CeRuSi 2 . 
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Introduction. Quantum phase transitions take place 
at zero temperature and are due to the zero point quan¬ 
tum huctuations. These fluctuations around the quan¬ 
tum critical point (QCP) display scale invariance both 
in space and in time, and result in the influence of the 
QCP over a finite range of temperature T. Hence the ef¬ 
fect of quantum criticality is detectable without actually 
reaching absolute zero T. Typically, some observables 
display diverging behaviour upon approaching the QCP. 

In itinerant electron systems, the criticality parameter, 
which determines the proximity to the respective QCP, 
may depend on the electron filling factor, the pressure, or 
the magnetic field (which is related to filling, through the 
Zeeman term). A related review can be found in ref. [Ij . 
Often these systems display non Fermi liquid behaviour, 
in the sense that e.g. the temperature dependence of var¬ 
ious quantities measured differs from the standard Fermi 
liquid (FL) one [l|. 

Our work is motivated b y a number of experiments on 
eight different materials |2l4l2l| , which display typical FL 
behaviour for appropriately low T. That is, quadratic 
in T resistivity and linear in T specific heat. However, 
the prefactors of these quantities appear to diverge as 
the respective QCP’s are approached. We show, via 
analytic diagrammatic calculations, that these facts can 
be consistently understood as arising from the exchange 
of relevant fluctuations among the quasi-particles. Our 
approach assumes that we deal with weakly disordered 
metallic systems. 


The model. We consider the Green’s function 


G^’^{k,e) 


1 

e - Ck + ep ± il‘^To 


( 1 ) 


with Cfe the quasiparticle dispersion, ep the Fermi energy 
and To the momentum relaxation time due to impurities. 
In the weak disorder regime [Hill epTo >1. Tq ^ is 
important as a regulator in our calculations. In fact, the 


characteristic FL dependence of Im E in eqs. m, 
(ITTl) is due to the finite ^. 

The dominant electron-electron interaction is assumed 
to be the fluctuation potential (or fluctuation propaga¬ 
tor) 


V{q,uj) = 


-iu}/[Dq^ + r) + q^Y -b i 


( 2 ) 


with g the coupling constant, ^ the correlation length 
and a measuring the distance from the QCP. The crit¬ 
icality parameter a depends on e.g. the magnetic field 
H, as in the systems of interest mentioned below, like 
a = /i®, h — \H/Ho — 1|, s > 0, where Ho is the critical 
field. The combination (Dq^ + r) in V{q,uj) is considered 
here for the first time. The factor Dq^ indicates disorder 
induced diffusion of the quasiparticles, with diffusion co¬ 
efficient D [l3, IHj . This factor could also originate from 
antiferromagnetic damping, but this origin would not be 
consistent with our treatment of ferromagnetic fluctua¬ 
tions - c.f. below. The factor r may originate from the 
inadvertent presence of elastic magnetic impurities in the 
samples [2l| , and should be equal to the relevant scatter¬ 
ing rate Tg^. It can also originate from the fermiology of 
a clean system (without disorder). qo=0 corresponds to 
ferromagnetic (FM) fluctuations, while finite qq to anti¬ 
ferromagnetic (AFM) fluctuations. 

In general, ^ and a are expected to be related through 
an equation of the type P, [l^ = dq a + fi T’", with 

do,f3,w constants. Below we will assume the Gaussian 
regime 


^ ‘^ = dq a , 


(3) 


with (3 = 0. However, for the purpose of our calculations, 
we will treat ^ and a as independent parameters, and 
consider the Gaussian regime relation in the final results. 
This procedure is entirely consistent, as can be seen from 
the details of the calculations below. 
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(a) 


with Kq = {Dq^ + r) (^^(g — qoY + ®)- We consider the 
limit e —>■ 0. The function A{x) = coth(a;) + tanh(—x) 
satishes A(x —>■ 0) = 1/x + 0{x) and A(|x| » 1) —?> 0. 
Hence we evaluate the integral over oj as 



(b) 

FIG. 1: Diagrams (a) for the self-energy and (b) for the free 
energy. The continuous lines are the fermion propagators, i.e. 
the Green’s function of eq. 0, and the dashed line is the 
fluctuation mediated interaction V{q,u)). 


Calculation of the scattering rate. We calculate 
the quasi-particle scattering rate as a function of T —>■ 0. 
For the self-energy S we use the relation [l^ (c.f. pg. 
183) 


/ OO 

duj Im G^{k — q,e — uj) 

-OO 

Im w) {coth(a;/2T) -|- tanh((e — w)/2T)} , (4) 


r2T 


duj 


-2T 


2T 


R{q,uj) 


OJ 


q 


8T^ 

lo 

q 


Rid) , (7) 


where R{q) = R{q,Lj ~ T), and 2T < Kq is implied. 


First we consider the case go = 0- We see that the 
dominant contribution to the remaining integral L over g 
comes from finite q > qi = \/a/f 


Ro f 1 

2B [a + C^g^ 




( 8 ) 


with qmax = 1/‘2 toVf, Ro = Rid) and B = aD - rf'^. 

We consider in detail two different limiting cases, 
namely r = 0, D > 0 and r > 0, = 0. 

For the case r = 0, D > 0 we have 


in order to calculate the scattering rate, which equals 
twice Im E - c.f. fig 1(a) for the corresponding Feynman 
diagram. 

We present the calculation of Im E(/c, e —>■ 0). In 3-D 
the integration over g is J2q — 1/(27’')^ / dq q^ J_j^dx, 
where x = cos 9, and 9 is the angle between the vectors 
k and q. From the factor Im G{k — q, e — oj), we have the 
integral over x 


T=^H] {dmax\ If 1 _ 1 \ 

gi J 2a\edlax + a edf+aj 

(9) 

This yields L oc 1/a^. Then, for the case r > 0, D = 0 
we have 

^ = ~^r{edlax+a~eWT^} ■ 


R{q,oj) = / dx Im G(fc — g, e — Cli) = / — 

J-i J-i si 


dx 


-b(c-b5x)2 Here, L ocl/(C2 ^). 


1 

b 


arctan 



— arctan 



where So = l/2ro, c = e — oj — e^ — q^/2m -I- cp and b = 
kq/m. To be specific, we assume a parabolic dispersion 
relation Ck = kf /2m. However, the precise form of the 
dispersion is not of particular importance. We have 

Im V{q,oj)=guj , (6) 

oj-^ -I- 


( 5 ) To estimate Ro we consider the relation arctan(x) — 
arctan(j/) = arctan((x—?/)/(l-|-x?/)), which yields i?(g) ~ 
arctan(2fcg/mSo)- The most relevant momenta are k ~ 
kp, thus yielding k/m = vp. Then we take g = qmax/‘2 
and we obtain Rq ~ arctan(l) = 7r/4. 

Next, we turn to the case of finite go and, as above, we 
consider g > gi = -y/a/^ and the limit a —>■ 0. 

We have for all r, D 


J 


do Rq 

2(Dgg -b r)^ a 


f 1 

< —pz arctan 

I Va 


L = 

(^) 



dq 


d Rid) 

(Dq^+r){a + eid-do)^)^ 


C(g- go) 
a + f^iq-qo)'^ 


+ o{VK) 


Qmax 

Qi 


( 11 ) 
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We now turn to 2-D. Again, 0 is the angle between the 
vectors k and q. From the factor Im G{k — q,e — uj), 
we have the angular integral l{q,u:) = sl+{^+bcosey 
with So,c,h as above. For small uj and k kp, e ^ 0 
we can make the approximation I{q,uj) ^ = Iq. As 

above, the integral over uj is duj (2T /uj) {uj Iq)/ + 

AT^) ~ 8T^ Iq/K'^. Then, we have the remaining integral 
L 2 over q in 2-D 

pQmax J 

" 4 W+r-)(a+V(g-9o)")2 ■ 

We notice that the integrand is very similar to the 3- 
D case of eq. (ED, and differing only in the factor Iq. 
Hence, we obtain the same scaling of Im E with a as in 
3-D. 

In 3 dimensions the final result is 

Im E(fcF,e = 0,T,a) = g/(27ruF)/3(a,go) , (13) 

with VF the Fermi velocity. The function / 3 (a, go) is given 
in Table 1. Therein Iq = qo/{DqQ + r). 

In 2 dimensions we have 

Im E(fcF,e = 0, T, a) = (8g To/f)/ 2 (o, go) . (14) 

It turns out that / 2 (a,go) = / 3 (a,go), so the prefactor 
dependence on a is the same as in 3-D. We note that the 
above dependence on a, ^ and T is valid for all k of the 
order of qmax or greater. For k away from kp only the 
prefactors change. 

A non-FL result Im E oc '/T was obtained in for 
dean metallic systems close to a FM QCP. We emphasize 
that herein we treat weakly disordered systems instead - 
c.f. the comment on following eq. (HD above. 


Calculation of the conductivity. We consider the 
total quasi-particle scattering rate 

r-i(r,a) =r-42ImE(fcF,e = 0,r,a) = 25' , (15) 


with T~l due to impurity scattering. 

We calculate the total conductivity a from an infinite 
series of diagrams involving disorder and V{q,uj). The 
n-th term of the conductivity series, shown in fig. 2, 
comprises n impurity scattering lines in parallel. We re¬ 
call that two impurity lines crossing each other introduce 
a small factor l/epTo -Cl - c.f. refs. MM . Hence we 
ommit any such diagram. In ref. we summed up to 
infinite order another diagrammatic conductivity series, 
which includes disorder and interactions (and yields ex¬ 
perimentally observed positive giant magnetoresistance). 
The Green’s function here is taken as 


Gf’^(fc,e) 


1 

e — Ck + ep ± iS ’ 


(16) 


i.e. it includes the self-energy of eq. m due to the 
fluctuation potential D(g, w). Due to momentum con¬ 
servation, the momenta in the upper and lower lines of 



FIG. 2: Feynman diagram for the n-th term of the conduc¬ 
tivity series. The dashed lines with crosses stand for impurity 
scattering. There are n such impurity lines inside the con¬ 
ductivity bubble. The Green’s function here is given by eq. 

(ED- 


the m-th pair of are the same. The major con¬ 

tribution to these diagrams comes from assuming for the 
various vectors ki = k 2 = k^ = ... = kn = kn+i. Also 
e = 0. 

Ommitting vertex corrections, the n-th term of this 
series is given by (cr = X)4o 

f dk e {Gf{k, 0)Gf{k, 0))"+i (17) 

J 

=u" ^ r dx 

° TO (CC2 + 52)-+l ’ 

where x = Ck — ep, Np is the density of states at the 
Fermi level, e is the charge of the electron, m is the mass 
of the electron, Uo = nimpV//^p, with nimp the concentra¬ 
tion of impurities, Vimp the typical value of the impurity 
scattering potential. Summing up this series, and con¬ 
sidering the energy scale Eq > ep, the result is 


‘ie'^Np 


ep 




arctan 


cf 




In 


4 + 

^ + 52 - ? 


.(18) 


For Eq ^ 0{eF) and ep > S we may approximate a = 
4 f e^Npcp/lm •\/52 — Uo). 

As the resistivity p = 1/a comes out proportional to 
the scattering rate, these results compare very favorably 
with data on CeCoIns Sr 3 Ru 2 07 Q and YbRh 2 Si 2 
Si , where p displays typical FL dependence 

p(r, a) = po + A(a) , (19) 


with po=const., and A{a) diverging around the QCP, i.e. 
around a = 0, as in our result. With all the results 
of Table 1 in mind, we may concretely assume that the 
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Gaussian 3-D case go = 0, r = 0, D > 0 applies (with the 
alternative being the 3-D case go = 0, r > 0, D = 0, and 
^ not scaling with a, as mentioned below). Assuming the 
scaling A{a) oc a“^, as in the case go = 0 (FM) above, 
the diverging A{a —>■ 0) in CeCoIns yields 2s = 1.29±0.1 
for the resistivity and 2s = 1.34 ± 0.1 for the thermal 
resistivity [1, Q. Hence, for CeColns, the Wiedemann- 
Franz law, i.e. a constant ratio of the thermal conductiv¬ 
ity K over the electrical conductivity cr, times T, in the 
low T limit, i.e. k,/ aT=coTist., is obeyed Q. These ex¬ 
perimental results are easily understood in the frame of 
our calculation: conduction electrons carry both charge 
and heat, while interacting via V (g, ui) (only small energy 
transfer is involved with this V{q,uj)). 

The same scaling A(a) oc was also found to fit 
data in La 2 -a:Cea;Cu 04 [Tl] with 2s = 0.38, in overdoped 
Tl 2 Ba 2 Cu 06 +a; @] with 2s = 0.62, in CeAuSb 2 @ with 
2s = 1.0, and in YbAlB 4 [l^ with 2s = 0.50. 

In principle, it should be possible to probe the scatter¬ 
ing rate through angle-resolved photoemission (ARPES) 
experiments (which require, of course, that the materials 
in question be adequately cleavable). 


Here the Matsubara energies are ei = {21 + 1 )ttT and 
LOm = 2rmrT, Dr, = v\Told (d=2,3 according to dimen¬ 
sionality) and [ij n(g,*a;m) = TJ2k,ei^i^ + + 

iuijn) G{k,iei) = Nf{ 1- To{Doq^ A l^^ml))- 

The integrand P(g, uj) = V (g, w) n(g, w) has a branch 
cut for Im w = 0, which comes from Im V (g, w) (in 
which —iuj —>■ -|-|a;m|)- We ommit the g, w dependence of 
n(g, w). Then, via Cauchy’s residue theorem, we obtain 
^ ImP(g,a;). 

The Bose distribution function nB{aj) = l/{e^^/'^'> — 1) 
is commonly approximated [2lj as nB{uj) = {T/lo) 6{T — 
|a;|). We first carry out the lo integration and then the 
integration over the momentum g in 2 and 3 D. 

Setting Yq = ^^(g — go)^ -|- a we obtain 



\UJ) {Dq^ + r) A /{Dq^ A rY 
2T { T \ 

Yq \{Dq^Ar)Yq) '' ’ 


Calculation of the specific heat. We calculate the 
specific heat C{T) = —Td‘^F{T)/dT‘^ = 7 T, through 
the free energy F{T) - c.f. fig 1(b) for the relevant Feyn¬ 
man diagram. To lowest order in V (g, w) we have 

F{T) = T‘^ ^ V{q,ii^m) G{kAq,KiAiuJm) G{k,iei) 

k,q,ei ,ujm 

= T ^ V{q,iuJm) U{q,iuJm) 

q,U)m 


We consider the low T regime T < {Dq^ A r) Yq and, 
as a result, F{T) satisfies 

F{T) = -r(T) . (22) 

First, we consider the case go = 0. In 3-D we have, with 
and qa = min | y^T/(aD), (r/(D^2))i/4|^ 


C3 = 

.(20) remaining integral over g 


r(T) = C3 


CQmax 


'Qa 


dq g2 


C3 


qB 


{Dq^ A r)Y^ 2B^ a -f 


— 2'/rD arctan 


'qVD\ 


——=— arctan —= 

CVa Vva 


r 




To proceed, we take the limit T —>■ 0 first and then the 
limit a —> 0 , in order the extract the coefficient 7 . 


Next we turn to the case of finite go. We set Ao = 
Dgg -I- r. Further, the minimum g = qmin should now 


Likewise, in 2-D, with C 2 = g we have for go = 0 satisfy qmin = \JTjD /(^ max{go, qmax})- In 3-D 


r(r) = C 2 


dq q 


Iq^ {Dq^ A r)Y^ 


(24) 


C 2 \ 1 1 . / Dq^ A r 

2B\aA Sfiq^ ^ B \a A Sfiq^ 


r(T) = C3 


dq g 2 




' Qmi 


{Dq^ A r)Yq^ 2^Xoa\^ 


V 

I- 


C(9 - qo) 

^Aeiq-qoY 


A O {Va) 


(25) 


And similarly for 2-D. 


The coefficient 7 in 3-D and 2-D is shown in Tables 1 
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and 2 respectively. Due to phase space considerations, 
for go = 0 we have 72 D oc 

These results are consistent with CeCoIns data Q 
as a function of a{H). They are also consistent with 
YbRh 2 Si 2 0 and Ge-doped YbRl^Si 2 data 0 (c.f. fig. 
2 therein), with Sr 3 Ru 207 data [ll|, and with CeRuSi 2 

data fI3l. 


Kadowaki-Woods ratio. The scaling of the 
Kadowaki-Woods (KW) ratio A/ 7 ^ in 3-D and 2-D is 
shown in Tables 1 and 2 respectively. Upon assuming the 
Gaussian regime, c.f. eq. ®, the KW ratio is constant 
only for the 3-D case go = 0, r = 0, D > 0, modulo the 
logarithmic in a divergence. However, if, alternatively, 
^ and a were independent parameters, and ^ would not 
scale with a, the KW ratio would be constant only for 
the 3-D case go = 0, r > 0, D = 0. In all other cases 
either the ratio goes to zero for small a, or there are no 
diverging prefactors H, 7 , contrary to the experiments. 

Also we mention that substituting the factor (Dq^ + r) 
in eq. ( 2 ) by (vo g), i.e. the usual ferromagnetic damping 
with Vo a constant, does not yield a (quasi-)constant KW 
ratio both in 3-D and 2-D. 

An a-independent KW ratio was observed in CeGoIns 
YbRh 2 Si 2 0, Ge-doped YbRh 2 Si 2 [1], and YbAlB 4 
10( 1 though (in most cases) in a more restricted range of 
H than the scaling of the coefficient A{a). E.g. in @ 
the KW ratio increases as the control parameter a —?> 0, 
which is consistent with the logarithmic dependence on 
a. 


For the other materials mentioned above, the experi¬ 
mental data are incompletely known with respect to the 
KW ratio. Therefore, the possibility exists that they fall 
in some other case, among the ones mentioned in Tables 
1 and 2, such that the KW ratio is not constant. 


In ref. [17|, using a different approach, diverging FL 
prefactors were obtained both for the resistivity and the 
specific heat. However, the results differ from ours. 
Therein, the KW ratio is constant only for 2-D FM fluc¬ 
tuations. 


A number of experiments, probing quantities other 
than the above mentioned, suggest AFM behavior in 
GeCoIns [ 2 , 23|. A possible explanation is that both FM 
and AFM fluctuations coexist in this material, with A{a) 
and 7 (a) being determined dominantly by FM fluctua¬ 
tions. 

We also consider an interaction with peaks at 
specific wavevectors goi (AFM case) K(g,w) = 
Er=i 5 /{-*w/(Dg 2 -hr)-h^ 2 (g-go*)^-ha}. In 2-D for 
tetragonal symmetry = 4 and in 3-D for cubic sym¬ 
metry nd = 6 . For small a, in 2-D and 3-D the potential 
U,(g,w) gives the same scaling of the prefactors as for 
the finite go cases above. 


Overview. In all, a consistent Fermi liquid descrip¬ 
tion emerges from these calculations. The renormaliza¬ 


tion of the fermions due to V (g, w) leads back to the FL 
fixed point in a low-T part of the phase diagram. This 
is consistent with experiments - e.g. c.f. fig. 3 of ref. 0 
for the case T^. The prefactors for the scattering rate, 
the resistivity and the specific heat diverge as power laws 
of the criticality parameter a. According to our calcula¬ 
tions, the Kadowaki-Woods ratio is constant only for 3-D 
FM (go = 0) fluctuations (possibly times the logarithmic 
in a divergence). 

We did not calculate explicitly the effective mass m* 
of the electrons. Experiments in Sr 3 Ru 207 [131 have 
indicated the absence of a magnetic field H dependent 
renormalization of to* (definitely so for five out of the six 
bands). This result is not inconsistent with our calcula¬ 
tions, where the diverging overall prefactors, as a func¬ 
tion of a, come from the small g dependence of the po¬ 
tential V{q,u}), and yield a uniform a dependence within 
the Fermi surface (c.f. the comment after eq. ca). 

Finally, we comment on the linear in T resistivity p 
displayed by Sr 3 Ru 207 [3] and CeGoIns in the vicinity 
of the QCP and for not very low T. In ref. [ 2 ^ we 
developed a fully microscopic FL model with a strong 
van Hove singularity (or, strong peak in the density of 
states), located at a characteristic energy e^H close to 
the Fermi level p,. This yields a quasi-particle scattering 
rate which is linear in T for T > (p — The 

model works very well for the cuprates, and this is how 
the linear in T resistivity of La 2 -a;Cea:Gu 04 Q can be 
understood, given the small difference p — CyH for many 
cuprates [26|, as shown by ARPES expts. Such expts. 
[271 ] on Sr 3 Ru 207 indeed yielded p — e^n = 4 meV, in 
agreement with our model [ 2 ^ . It is fair to attribute the 
linear in T resistivity of GeCoIns to the same mechanism, 
i.e. originating from a van Hove singularity, which resides 
close to the Fermi surface in the vicinity of the QGP. 
Moreover, it is reasonable to view the regime p oc T^, 
I < 6 < 2, displayed by Sr 3 Ru 207 [3] and GeCoIns [^, 
as a smooth transient between 6=1 and 6 = 2 . 

* e-mail : kast@iesl.forth.gr ; giwkast@gmail.com 
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ae 
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Xq a \qmax-qo qmin-qo J 

Qmax ^ QO ^ Ql ^ 0 

In TT 

2 £ a3/2 

=3 9n 

__ 
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Supplementary Information 

Appendix A : On the effective interaction 

V{q,U}) 

V{q, oj) of eq. ( 2 ) can only be derived in the context of 
an RPA-type approach, in the spirit of references [15-20] 
cited in the article. This fact is also emphasized in the 
recent article by Y. Wang and A.V. Chubukov in Phys. 
Rev. B 92, 125108 (2015). 


Appendix B : On the calculation of the scatter¬ 
ing rate 

In the limit T — ?> 0 the thermal function X = 
coth(w/2T) -I- tanh((e — u})/2T) in eq. (4) becomes 
A = 2 for 2T < w < e, and A = 0 for w < —2T and 
uj > e. Then the integration over w - compare with eq. 
(7) - amounts to 


~ g i?o In 


2 


/ du! Im V (g, w) R{q, oj) 

J 2 T 





(26) 


for Kq > e. The rest of the algebra proceeds as in eq. 
( 8 ) and onwards. Thus the scattering rate scales like 
as well, as expected for the FL regime. 


We also give some details of the 2-D calculation in the 
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main text. The angular integral I{q,uj) is given by 

d9 s„ 


= [ 
Jo 


0 + (c + 6cos6<)2 


1 


+ C.C. , (27) 


y 1/52 — c? — 2 iSoC + Sq 
with c, b as above. Then, taking 
a = + si > 0, , (3 = 2soC , a » |/3| , (28) 

we can approximate 


I{q,uj) ~ 27r 


, 1/2 o 

a \ ZTT 


a2 + /32 


(29) 


Appendix C : On the calculation of the infinite 
series for the conductivity 

Setting Zo = ANpe^/m, the n-th term of the conduc¬ 
tivity series is given by 


A„ = m” Zo / dx- 


ep + X 


(x2+^2)"+l 
= Mg Zo {ep r„ — ep An-\- Bn) , (30) 


where x = Ck — ep and 

dx TT (2n — 1)!! 


r„ = 


(^2 S' 2 )n +1 2" S'2"+1 n! 


(31) 


/-oo (x2+^2)n+l ’(32) 


Bn = 


dx : 


(x2 + ^2)n+l 


.(33) 


We easily obtain 

00 

r = ^< r„ = 


n—0 


-Uo 


(34) 


A=J2< ^n = 


1 


n—0 


- Uo 


— — arctan ( , ^ = 

2 WS'^-Uo 


(35) 

For n = 0 only^ we consider as the upper limit of inte¬ 
gration the energy scale Eq > ep, instead of infinity (this 
integral is ultra-violet divergent), and we obtain 


^°=2 


1 dy 


y + S^ 


1 fEl + S^ 


(36) 


For n > 1 we take infinity as the upper limit of integra¬ 
tion, thus obtaining 


Bn = bn 


Hence 


B = Y,u:Bn = -\n 


n—0 


Cp -I- 52 - Uo 


El+ 8'^ 

Cp -I- 52 - Uo 


(37) 


(38) 


Putting all these terms together yields the total con¬ 
ductivity given by eq. m above. 
























